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It  is  shown  that  several  families  of  PBIB  designs  are 
E-optimal  over  the  collection  of  all  block  designs.  Among 
these,  the  partial  geometries  with  two  associate  classes^/u: 
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and  either  i  -  1  <  Jv  <  k  or  k  <  /v  <  i  -  1)  r^The  duals  of 

these  designs  are  also  E-optimal.  When  uniform  on  rows,  these 
designs  remain  E-optimal  in  the  additive  setting  of  two  way 
elimination  of  heterogeneity. 


Key  words  and  phrases:  E-optimality,  partial  geometry,  PBIB 
design,  information  matrix,  efficiency,  eigenvalue. 

AMS  1970  subject  class i f ica tion :  Primary  62K05,  Secondary  62K10. 

Research  silpporte^d  in  part  by  the  Air  Force  Office  of  Scientific 
Research  under  the  grant  AFOSR  80-0170 


1 


1.  Introduction  and  notation 

The  object  of  this  work  is  to  investigate  the  E-opti¬ 
mality  of  discrete  statistical  experiments  in  the  additive 
setting  of  one  and  two  way  elimination  of  heterogeneity.  For 

v  ,  b  ,  k  positive  integers,  we  denote  by  Q  the  col- 

v,o,k 

lection  of  all  k  x  b  arrays  with  varieties  1,  2,...,v  as 
entries  (2  <  k  <  v)  .  Any  such  array  d  e  ft  ,,  is  called 

V  ,  D  ,  K 

a  design.  The  columns  of  d  are  called  blocks.  A  design  is 
said  to  be  binary  if  each  block  of  d  consists  of  distinct 
varieties;  d  is  called  equireplicated  if  each  variety  occurs 
the  same  numbers  of  times  throughout  the  whole  array  d  .  If 
each  variety  appears  the  same  number  of  times  in  each  row  of 


d  ,  d  is  called  uniform. 

Let  be  the  unknown  effect  of  variety  i  and 

be  the  (unknown)  effect  of  the  j ^  block.  In  the  additive 

model  of  elimination  of  heterogeneity  in  one  direction  /  we 

assume  that  the  expectation  of  an  observation  on  variety  i 

in  the  block  of  d  is  a.  +  ft  -  • 

J  x  H  j 

If,  in  addition,  Y^  denotes  the  (unknown)  effect  of  the 
m^^  row  of  d  ,  and  there  is  evidence  of  row  effects  in  the 
observations,  we  shall  assume  that  the  expected  value  of  an 
observation  on  variety  i  in  block  j  and  row  m  of  d  is 

ai  +  6j  +  Ym  * 

Under  both  models  we  assume  the  kb  observations  uncor- 

2 

related,  with  common  (unknown)  variance  o  .  The  main  in¬ 
terest  is  in  comparing  the  variety  effects  • • •  *av  * 

The  information  matrices  of  variety  effects  under  the 
two  models  are,  respectively 


and 


kCd  k  tli  ay  (rfU  , .  .  .  ,  rdv)  -  Nj 


kCd  =  k  diag(rdl . r^)  -  -  b“  Md<kI-J)M£ 

where  Nd  •=  (ndij)  ,  Md  =  (m^)  .  with  ndij  (resp.  lUj..) 

indicating  the  number  of  times  variety  i  appears  in  the 


j  block  (resp.  row)  of  d  ;  is  the  replication  number 

of  variety  i  in  d  .  J  denotes  the  matrix  with  all  its 

entries  1  and  I  is  the  identity  matrix.  By  we  denote 

the  (i/j)*"*1  entry  of  *  It  is  known  that  for  any  d 

both  and  are  nonnegative  definite  with  row  sums 

zero.  Let  further  0  =  y  ,A  <  u,  ,  and 

^dO  dl  Kd#v-1 

0  =  {J  j  n  <  p,,  <...<  ij  ,  ..  be  the  eigenvalues  of  C,  and 

dOdl  d ,  v- 1  ^  d 

,  respectively. 

A  design  d*  is  called  E-optimal  over  ft  ,  ,  (under 

V  |  D  ,  K 

a  given  model)  if  the  maximal  variance  of  normalized  best 
linear  unbiased  estimators  of  estimable  functions  is  minimal 
under  d*  .  It  is  well-known  (Ehrenfeld  (1955) ;  see  also 
Kiefer  (1959)  and  the  beautiful  expository  article  of  Kiefer 
(1978))  that  d*  is  E-optimal  over  ft  ,  ,  under  the  one  way 

V  ,  D  ,  X 

elimination  (resp.  two  way  elimination)  if  and  only  if 

yd*l  -  wdl  (resp-  yd*l  -  ydl>  for  311  d  e  “v,b,k  *  In 
short.  E-optimality  deals  with  the  association 

d  -  cd  -  v,di 


and  the  objective  of  finding  the  design  d  witn  maximal 

d  J 

The  experimental  setting  which  involves  relatively  few 
blocks  (say  0  <  b  <  2v)  is  of  notable  practical  importance. 
Just  to  focus  attention,  suppose  we  have  v  -  15  varieties. 
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b  s  15  blocks,  k  -  3  varied  ips  per  block  and  that  there  is 
evidence  of  an  additive  variety  and  block  response*  The  ex¬ 
perimenter  would  like  to  use  a  design  d°  which  (if  not 
optimal)  is  efficient  in  the  E-scnso ,  i«e. ,  ^d°l  ^  ^dl 
most  d  .  You  would  then  say,  well,  maybe  this  d°  will  do: 


1  1  8  8  12  14  1  13  13  9  9  11  2  13  15 

d°  :  14  10  10  4  10  4  6  2  7  7  3  7  3  5  12  (1. 

8  4  6  12  14  6  15  9  3  5  11  2  5  15  11 

It's  connected,  it's  binary,  it’s  equireplicated ,  and  a  pair 
of  distinct  varieties  appears’  in  at  most  one  block.  Looks 
good,  as  it  is  M  -  S  optimal,  (see  Eccleston  and  Hedayat  (1974) 
or  Shah  (I960))  and  surely  does  not  favor  any  variety  (in  the 
senses  mentioned  above).  Only  trouble  is,  d°  is  in  fact  a  very 
inefficient  design  in  the  E-sense. 

With  a  bit  of  luck  you  would  have  probably  recommended: 


1 

1 

1 

2 

2 

2 

5 

5 

5 

8 

8 

9 

9 

7 

6 

d* :  4 

10 

3 

4 

3 

10 

4 

3 

10 

7 

12 

6 

11 

11 

14  (1.2) 

7 

13 

6 

9 

8 

14 

12 

11 

15 

15 

13 

15 

13 

14 

12 

which  has  features  very  similar  to  d°  in  the  M  -  S  sense 
(but  luckily  not  in  the  E-sense) . 

However,  one  should  hope  that  you  haven't  sent  him  back 
with  such  a  design  as 


1 

1 

1 

1 

i 

1 

i 

1 

i 

1 

1 

1 

10 

11 

13 

d:  1 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

10 

11 

13 

15 

1 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

4 

12 

14 

8 

—  • 


A 


The  .surprizing  fact  is  l  hat  d  is  E-boLLor  than  d°  and 
the  discrepancy  in  the  H-perf orinance  between  d*  and  d°  is 
by  a  factor  of  more  Lhan  8.  The  exact  figures  are  3y^0^  “  *571  , 

3ydl  =  and  3yd*l  =  5  * 

We  would  like  to  stress  two  things.  That  this  is  not 

an  isolated  example  (not  when  the  number  of  blocks  is  small) ; 

and  that  M  -  S  optimal  designs  with  large  number  of  blocks 

do  not  differ  in  the  E-performance  in  such  a  surprizing  way. 

In  the  former  case  a  way  to  avoid  bad  M  -  S  optimal  designs 

is  suggested  in  Constantine  (1980). 

In  the  next  section  we  show  that  a  design  d*  e  Q,  ,  . 

v  /  d  ,  k 

with  1^*1  ~  v^k  (k-1)  and  r  =  bkv  ^  integral,  is 

E-optimal  over  all  block  designs.  This  result  has  a  number 
of  consequences  in  settings  with  b  small.  Various  families 
of  P BIB  designs  (mentioned  below)  with  X^  and  ^  zero  or 
one  are  proved  E-optimal  over  all  block  designs.  The  partial 
geometries  with  two  associate  classes  are  such  instances. 

Among  the  triangular  PBIB  designs  with  schemes  of  size  n: 
those  with  =  0  r  X2  =  and  those  w^th  X^  =  1  , 

X2  =  0  and  k  >  n  -  1  .  Among  the  ones  with  schemes: 

those  with  X^  =  0  ,  X^  =  1  ,  k  ^  /v  and  those  with  X^  =  1  , 

X2  =  0  and  either  i  -  1  <  /v  <  k  or  k  ^  /v  £  .i  -  1  •  The 
duals  of  these  designs  are  E-optimal  as  well.  These  E-optimality 
results  extend  easily  to  the  setting  of  two  way  elimination  of 
heterogenci ty  when  the  above  designs  are  also  uniform. 

Results 

We  denote  by  1  the  column  vector  with  all  its  entries  1 


2. 
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and  by  J  the  (not  neccssa r i  l.y  squire)  matrix  with  nil  its 

entries  1.  The  following  two  lemmas  provide  upper  bounds  for 

u  .  Various  bounds  can  also  be  found  in  Chakrabarti  (1963) 
Hdl 

Jacroux  (1980),  Cheng  (1980)  or  Constantine  (1979b). 

Lemma  2.1:  Let  C  be  a  v  x  v  nonnegative  definite  matrix 
with  zero  row  and  column  sums.  Denote  the  eigenvalues  of 

C  by  0  -  <  y ^  <...<  *  Then  the  sum  of  entries  in 

-  _  ,  ..  ,  m (v-m) 

any  m  x  m  principal  minor  of  C  is  at  least  — — - 

(1  £  m  £  v  -  1)  . 


Proof:  Observe  that  a  matrix*  obtained  form  C  by  row  and 

(same)  column  permutations  has  the  same  eigenvalues  as  C  . 
It  will  therefore  be  enough  to  prove  the  lemma  for  the  mxm 
leading  principal  minor  of  C  .  Call  this  leading  principal 
minor  M  .  Then  . 


is  stated.  The  inequality  relies  on  the  known  fact  that 


V 


1 


=  min 
x  1 1  =  0 


x  [  Cx 
x  1  x 


and  on  observing  that  ^  Qj J  -  1^  1  =  0  (since  the  1  in 

is  m  x  1)  .  This  ends  the  proof. 

Our  next  lemma  gives  an  upper  bound  for  y^  when 
d  c  Q  ,  K  v  is  equireplicated. 

Lemma  2.2:  If  an  equireplicated  design  d  e  k  contains 

a  block  which  consists  of  m  distinct  varieties  (2  s  m  £  k) 


then 


kp 


dl  m(v-m) 


(k-1) (mr-k)  . 


Proof :  By  eventually  relabeling  the  varieties  and  reshuffling 

the  blocks,  we  can  assume  that  the  first  block  in  d  consists 


°f  ndU  l's  ,  ndn  2 '  s 


and  n.  .  m's  .  Index  the 
dml 


rows  and  columns  of  c,  by  the  varieties  1,  2,...,v  (in 

d 

this  order) ,  and  let  be  the  m  x  m  leading  principal 


=  r  and 

Hence 


minor  of  C  ,  .  Observe,  firstly,  that  Y  n _ .  . 

d  jhi  dij 

b  2  t>  2  b  3_i 

that  7  n, .  .  2  n , . .  +  Y  n  , .  .  >  Y  n , .  .  =  r  . 

£  di]  dil  dii  dii 

3=1  3=2  J  3=1  J 

^2  ^2 
■  Lndij  *  r  -  ndil  and  therefore  .  Y  ndij  *  ndil  +  r  "  ndil  * 
J=z  1=1  J 

J  m  J 

Secondly,  note  that  Y  A  , .  .  (a  sum  of  m(m-l)  nonnegative 

i?j  013  . 

1  J  m  b  m 


terms ) 


m 

)  satisfies  Y  \  . .  =  Y  Y  n , .  n , .  >  Y  n , . ,  n...  . 

i^j  dl3  ih  u=i  dlu  d3U  iii  dl1  d31 


m 


Using  these  two  inequalities  and  the  fact  that  £  n  =  k  , 

i=l  dl1 

we  obtain: 

m  b 


1  Md  1  =  mrk  -  Y  Y  n, . .  -  )  \  a .  . 

.  *■,  .  dii  . 4 .  dii 
i=lj=l  J  l^j  J 


m 


m 


<  mrk 


y  (n. .  ,  +  r  -  n  , . ,  )  -  Y  n,  .  ,  n-,.. 
dil  dil7  dil  dll 


m 


m 


=  mrk  -  (  Y  ndU)  -  mr  +  [  n  . 

1=1  1=1 


=  (k-1) (mr  -  k) 

That  ky,.<  — 7 — — — r  (k-1)  (mr  -  k)  follows  now  from  Lemma  2.1. 
dl  m(v-m) 

This  ends  the  proof. 

Through  the  remainder  of  the  paper,  let  the  varieties  in 
a  design  d  e  ft  .  ,  be  always  labeled  so  that  the  replication 
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numbers  rdi  satisfy  rdl  <  rd2  <...<  rdv  .  We  are  now 
ready  to  prove  our  first  result. 

Theorem  2.1:  Let  r  =  bkv  ^  be  an  integer.  A  design 

d*  c  Qv  b  k  which  satisfies  k^]d*x  *  (r-1)  (k-1)  is 

E-optimal  over  all  block  designs. 


Proof:  Let  d  be  any  design  in  Q__  ^  .  Then  d  is  either 

~  V  ,  D  /  K 

equireplicated  or  it  is  not.  Suppose  it  is  not.  Then 
r^  <  r  -  1  and  by  Lemma  2.1  with  m  =  1  we  have: 


k  y 


dl 


v 

v-1 


V 


rdl(k_1)  5  v-1 


(r-1) (k-1)  < 


^(r-1)  (k-1)  s  kudVtl 


which  show  that  such  a  design  is  strictly  E-worse  than  d*  . 

Assume  now  that  d  is  equireplicated.  We  may  also  assume 
that  d  has  a  block  which  consists  of  m  distinct  varieties 
(2  <  m  <  k) .  (Observe  that  if  d  has  no  such  block f  the  infor¬ 
mation  matrix  Cd  is  then  the  zero  matrix,  and  hence  for  such 
d  we  have  y^  =  0  <  •)  By  Lemma  2.2  we  can  write 

kpdl  -  mjv-mj  (k'1}  {mr  ‘  k)  ' 

2 

Let  Q(m)  =  -ky^^^m  +  (vky^*^  -  v(k-l)r)m  +  vk(k-l)  .  Note 

that 


m(v-m)  (k~1}  (mr  "  k)  -  kvd*l  '  for  a11  2  <  m  <  k 


if  and  only  if 
is  a  quadratic 
Q ( 0 )  =  vk (k-1 ) 
that  Q(m)  >  0 


Q (m)  >  0  ,  for  all  2  <  m  <  k  .  Since  Q  (m) 

in  m  with  negative  leading  coefficient  and 
>  0  checking  that  Q(k)  >  0  would  insure 
for  all  2  <  m  <  k  .  By  assumption 


'•W-  - 


8 


2 

kPd*l  "ijj-(r-l)  (k-1)  which  implies  -k  +  vk''d*i  “ 

v(k-l)r  +  v (k-1 )  >  0  .  In  terms  of  Q  this  last  inequality 
is  simply  k_1Q(k)  s  0  .  Since  k  is  positive  it  follows 
that  Q(k)  >  0  ,  as  desired.  We  have  therefore  shown 

k"dl  k  <mr  '  k>  5  kgdn  '  for  a11  2  5  m  5  k  • 

This  concludes  the  proof. 

The  following  aspect  of  Theorem  2.1  should  perhaps  be 

stressed.  It  is  easy  to  show  that  for  two  (real)  matrices 

A  and  B  ,  the  products  A B  and  BA  have  the  same  set  of 

eigenvalues.  Given  an  equireplicated  design  d  e  ^  y 

(with  r  =  bkv  1  integral),  we*  call  6  e  St  the  design 

d,  v,r 

dual  to  d  ,  if  .  Hence  rC^  =  rkl  -  *  BY 

the  above  remark  it  directly  follows  that  kC,  and  rCr 

a  6 

have  the  same  set  of  eigenvalues.  Thus,  if  d*  ^  ft  ,  . 

V  ,  o  ,  Jc 

satisfies  the  assumption  of  Theorem  2.1,  then 


rp 


6*1 


=  kp 


d* 1  v-k 


(r-1 )  (k-1)  - 


b-r 


(k-1) (r-1) 


and  hence  6* e  ft.  satisfies  also  the  assumption  of  the 

d  ,  v  /  r 

theorem.  This  gives 


Corollary  2.1  Let 

r  =  bkv  ^  be  an  integer. 

The 

dual  of  a 

design  d*  €  ft  .  . 

*  v ,  b ,  k 

with  vd*1  >  v_k  (r  1)  (k-1) 

is 

E-optimal 

over  ft. 

b ,  v ,  r 

All  the  designs 

which  we  shall  prove  E-optimal 

next,  have 

E-optimal  duals,  in  view  of  this  corollary. 

Although  not  apparent  from  the  statement  of  Theorem  2.1, 
the  lower  bound  ^-(r-1)  (k-1)  is  useful  for  designs  with 
relatively  small  number  of  blocks  (0  <  b  <  2v  ,  say),  and 


hardly  over  otherwise, 
of  Theorem  2.1. 
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hot  u  e>:nini  no  a  fow  consequences 


Theorem  2 . 2  A  Partially  Balanced  Incomplete  Block  design 
based  on  v  varieties  and  b  blocks  of  size  k  ,  with  in¬ 
tegral  parameters  as  below,  is  E-optimal  over  all  block  de¬ 
signs:  r  =  bkv  1  ,  ^  =  1  ,  =  0  ,  nL  =  bk(k-l)v_1  , 

n2  =  v  -  1  -  bk (k-1 ) v-1  ,  t  =  k  (k-1)  (bk-v) v_1 (v-k)-1  , 

p1^  =  k  -  2  +  [k  (k-1)  (bk-v)  -  v  (v-k)  ]  (bk-v)  v  ^ (v-k)  1  , 

P11  =  bk2(k-1)  (bk-v) v"2 (v-k) ~1  . 

Proof :  Let  a  PBIB  design  with  parameters  as  above  be  denoted 

by  d*  .  Then  is  known  to  have  two  distinct  nonzero 

eigenvalues  0  <  2  •  whore  kU^*^  =  r(k“l)  -  k  +  t  +  1  . 

For  a  proof  see  Connor  and  Clatwortny  (1954),  Bose  and  Mesner 
(1959)  or  Raghavarao  (1971,  p.195) .  It  is  straightforward  to 
check  that 

r  (k-1 )  -  k  +  t  +  1  =  kyd*x  =  ~j^(r-l )  (k-1)  . 

V7e  are  then  done  by  Theorem  2.1. 

The  partial  geometries  (r,k,t)  with  two  associate 
classes,  defined  by  Bose  (1963) ,  satisfy  the  conditions  of 
Theorem  2.2.  Whence, 

Corollary  2.2.  Partial  geometries  with  two  associate  classes 
are  E-optimal  over  all  block  designs. 

For  a  thorough  reading  on  classical  finite  geometries  we 
refer  the  reader  to  the  book  of  Dembowski  (1968) .  Many  PBIB 
designs  can  be  constructed  from  partial  geometries.  The 
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varieties  would  be  the  points  and  the  blocks  would  be  the  lines 
of  the  geometry.  We  mention  the  work  of  Seiden  (1961)  and 
Ray-Chaudhuri  (1962)  in  connection  with  this.  The  design  (1.2) 
is  an  example  of  a  partial  geometry  of  the  simplectic  type. 

Bose  and  Clatworthy  (1955)  showed  that  PBIB  designs  with 
b  <  v  ,  X^  =  1  and  X^  =  0  necessarily  have  parameters  as 
those  listed  in  Theorem  2.2.  We  therefore  have: 

Corollary  2 . 3.  A  Partially  Balanced  Incomplete  Block  design 
with  b  <  v  ,  X  ^  =  1  and  X ^  =  0  is  E-optimal  over  all 
block  designs. 

Connor  and  Clatworthy  (1954)  found  the  nonzero  eigen¬ 
values  of  the  information  matrix  of  a  PBIB  design  with  two 
associate  classes  to  be 

kw,  =  r(k-l)  +  +  Ai  +  A2} 

and 

kp2  =  r(k-l)  +  j[(X1~X2)(-Y  -  /A)  +  +  X2)l  • 

It  is  easy  to  see  that  y ^  if  and  onlY  if  <  ^2  # 

Y  and  A  are  expressed  in  terms  of  the  parameters  of  the 

2  1  2  12 
association  scheme  as  Y  =  p^  anc^  ^  =  ^12  ~  ^12^  + 

2 (P^2  +  p12^  +  1  (See  Ra9havarao  (1971,  p.  126)). 

We  now  prove  the  following; 

Theorem  2 . 3  (a)  A  Partially  Balanced  Incomplete  Block  design 

with  A.=  0  ,  A  -  i  and  Y  -  /A  +  1  *  is  E-opti- 

mal  over  all  block  designs;  (b)  A  Partially  Balanced  Incomplete 
Block  Design  with  and  1  -  7  -  /K  *  —  — Jr._ 
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is  E-optimal  over  all  block  designs. 

Proof :  Let  d*  be  a  PBIB  design  as  in  (a) .  Then 

Y  -  /A  +  1  >  2  ^  ~  implies  (v-k)  (y  -  /A  +  1) 

V  K 

2 

>  2rk  -  2vk  -  2rk  +  2v  which,  m  turn ,  implies 

r  (k-l )  (v-k )  +  -^-(v-k)  (Y  -  /A  +  1)  >  v(r-l)(k-l) 

and  which  can  be  rewritten  as 

r  (k-l)  +  j(Y  -  /A+  1)  >  (r-l)  (k-l)  . 

We  are  now  done  by  Theorem  2.1,  since  it  follows  from  the 
paragraph  preceding  the  statement  of  Theorem  2.3  that 
kPd*^  =  r(k-l)  +  -  /A  +  1)  .  As  for  part  (b)  ,  it  can  be 

similarly  proved. 

When  specialized  to  various  known  association  schemes 
Theorem  2.3  yields  several  corollaries: 

Corollary  2.4:  A  Partially  Balanced  Incomplete  Block  design 
with  a  triangular  association  scheme  of  size  n  ,  =  0  , 

X2  =  1  and  block  size  k  >  j  ;>  3  is  E-optimal  over  all 

block  designs. 

Proof ;  In  this  case  y  =  n  -  5  and  A  =  (n-2)^  ,  v  =  n 

and  r (k-l)  =  (P~2) jn~3)  .  Now,  for  n  >  6  we  have  the  fol¬ 

lowing  chain  of  implications:  k  ^  i  =£>}c2  ---  ^2Vjj — -j— ==> 

2v  s  (v-r  (k-l)  +  l)k=^vk  -  v  >  r(k-l)k  -  v(k-l)  ==> 

Y  -  /A  +  1  =  -2  ^  V^~  »  as  desired.  The  proof  can 

now  be  concluded  by  part  (a)  of  Theorem  2.3. 

The  family  of  PBIB  designs  mentioned  by  Shrikhande  (1965) 
satisfies  the  assumptions  of  Corollary  2.4.  These  PBIB  designs 


_ 
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are  therefore  E-optimal  over  all  block  designs  (with  the  same 
v,  b  and  k)  . 

In  a  similar  way,  and  with  the  help  of  Theorem  2.3(b), 
we  can  prove 

Corollary  2 . 5 .  A  Partially  Balanced  Incomplete  Block  design 
with  a  triangular  association  scheme  of  size  n  ,  =  1  , 

=  0  and  block  size  k  >  n  -  1  is  E-optimal  over  all  block 
designs . 

Examples  of  such  families  can  be  found  in  Masuyama  (1965) 

Some  can  be  obtained,  for  example,  by  just  writing  the  rows  of 

the  triangular  association  scheme  as  the  blocks  of  the  design. 

We  now  turn  to  PBIB  designs  with  L.  schemes.  Here 
2 

v  =  s  (for  some  integer  s  >  2)  ,  Y  =  s  -  2i  +  1  and  A  = 
In  much  the  same  way  as  Corollary  2.4  was  obtained,  we  are 
led  to 

Corollary  2.6.  A  Partially  Balanced  Incomplete  Block  design 
with  an  association  scheme,  =  0  ,  =  1  and  block 

size  k  ^  /v  is  E-optimal  over  all  block  designs. 

When  /v  is  a  prime  or  a  prime  power,  Clatworthy  (1967) 
gives  a  class  of  PBIB  designs  with  the  following  parameters: 

2 

v  =  s  ,  b  =  s  (s  +  1  -  i)  ,  k  =  s,r  =  s  +  l-  i  ,*A^  =  0  and  A^  =  1 

By  our  corollary,  all  these  designs  are  E-optimal. 

Let  us  now  derive: 

Corollary  2.7.  A  Partially  Balanced  Incomplete  Block  design 
with  an  association  scheme,  A^  =  1  ,  A ^  “ 


0  and  the 
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block  size  k  satisfying  either  i  -  1  <  /v  <  k  or  k  <  /v  <  i  - 
is  E-optimal  over  all  block  designs. 

Proof :  With  Y  and  A  as  in  the  paragraph  preceding  Corol¬ 
lary  2.(7  we  have  -  y  -  SE  +  1  =  2  (i-s)  and  r(k-l)  =  is  -  i  . 

The  condition  i  -  1  <  s  <  k  or  k  <  s  <  i  -  1  can  be  written 
as  (s  -  i  +  1)  (k  -  s)  >  0  .  Upon  multiplying  by  s  and  ex¬ 
pending,  it  leads  to  vi  -  vs  +  ks  +  vk  -  v  -  kis  ^  0  .  Using 

the  fact  that  r(k-l)  =  is  -  i  this  gives  (v-k)  (i-s)  > 

(k-1)  (rk-v)  or 

-  y  -  /A  +  1  =  2  (i-s)  >  2-i?iri2  (rk-v) 

v— k 

With  this  last  condition  satisfied,  we  are  done  by  Theorem 
2.3(b). 

As  an  example  of  PBIB  designs  whose  parameters  satisfy 
the  assumptions  of  Corollary  2.7  we  mention  the  ones  with 

2 

v  =  s  ,  b  =  is,k  =  s,r  =  i,  =  1  and  A2  =  0  ?  (s  >  2)  . 

These  designs  appear  in  Clatworthy  (1967). 

Let  us  now  turn  our  attention  to  the  setting  of  two  way 

elimination.  Note  that  the  matrix  C ,  -  C -  «  k~*b~*M, (kl  -  J)M* 

ad  a  d 

is  nonnegative  definite.  Hence  (see  Bellman  (1979))  ydi  s  ydi  , 

1  s  i  s  v  *  1  .  If  d*  is  a  uniform  design  ,  then  M.  =  bv~^J 

d  * 

and  it  is  easily  seen  that  =  C  .  Moreover,  if  such  a  uni- 

d*  d* 

form  design  d*  is  E-optimal  in  the  setting  of  one  way  elimi¬ 
nation,  then  for  any  other  design  d  we  have 

udl  *  Mdl  S  wd*l  =  **d*l 

and  hence  d*  is  also  E-optimal  in  the  setting  of  two  way 
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elimination.  This  simple  argument  extends  to  D  -  ,  A  -  , 

type  I  criteria  (Cheng  (1978))  and  more  generally  to  Schur- 

optimality  (Constantine  (1979a).  Our  observation  can  be  for¬ 
mulated  as 

Lemma  2.3 :  A  uniform  E-optimal  design  in  the  setting  of  one 
way  elimination  is  also  E-optimal  in  the  setting  of  two  way 
elimination . 

The  results  obtained  so  far  do  then  extend  in  the  fol¬ 
lowing  way: 

Corollary  2.8:  if  any  of  the  designs  in  our  previous  theorems 
and  corollaries  are  uniform,  they  are  then  E-optimal  for  the 
settings  of  both  one  and  two  way  elimination  of  heterogeneity. 


It  is  often  possible  to  rearrange  the  varieties  in  each 
individual  block  and  thus  achieve  uniformity  of  the  design. 
As  an  example,  the  partial  geometry  (1.2)  can  be  written  as 
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d**:  3 
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Let  us  remind  the  reader  that  if  there  is  no  evidence  of 

row  effects  in  the  observations,  then  d*  and  d**  are  equally 

efficient  to  use.  But  if  row  effects  are  truly  present,  the 

E-performances  of  d*  and  d**  are,  respectively,  +  , 

d*l 

and  =  5  .  Trusting  that  the  response  is  indeed  additive, 

d**  is  then  clearly  to  be  preferred. 
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